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THE COMPUTATIONAL HARDNESS OF COUNTING 
IN TWO-SPIN MODELS ON rf-REGULAR GRAPHS 

* ALLAN SLY AND tNIKE SUN 



Abstract. The class of two-spin systems contains several important models, including 

random independent sets and the Ising model of statistical physics. We show that for 

both the hard-core (independent set) model and the anti-ferromagnetic Ising model with 

jrt ■ arbitrary external field, it is NP-hard to approximate the partition function or approximately 

sample from the model on d-regular graphs when the model has non-uniqueness on the 

d-regular tree. Together with results of Jerrum-Sinclair, Weitz, and Sinclair-Srivastava- 

Cs| ' Thurley giving fpras's for all other two-spin systems except at the uniqueness threshold, 

this gives an almost complete classification of the computational complexity of two-spin 

systems on bounded-degree graphs. 

^^ , Our proof establishes that the normalized log-partition function of any two-spin system 

n , ' on bipartite locally tree-like graphs converges to a limiting "free energy density" which 

coincides with the (non-rigorous) Bethe prediction of statistical physics. We use this result 

to characterize the local structure of two-spin systems on locally tree-like bipartite expander 

3 . graphs, which then become the basic gadgets in a randomized reduction to approximate 

MAX-CUT. Our approach is novel in that it makes no use of the second moment method 

employed in previous works on these questions. 

^ ■ 1. Introduction 

(N 

^ . Spin systems are stochastic models defined by local interactions on networks. The class of 

^ I spin systems includes well-known combinatorial counting and constraint satisfaction prob- 

^ ■ lems. In this paper we classify the complexity of approximating the partition function for 

O ' all homogeneous two-spin systems on bounded-degree graphs. 

^ ' When interactions favor agreement of adjacent spins, the model is said to be ferromagnetic. 

Jerrum and Sinclair [11] gave a fully polynomial-time randomized approximation scheme 
(fpras) for approximating the partition function (the normalizing constant in the probabil- 
r> ■ ity distribution) of the ferromagnetic Ising model, which covers all ferromagnetic two-spin 

c^ ! systems. For anti-ferromagnetic systems such as the hard-core and anti-ferromagnetic Ising 

models, the complexity of approximating the partition function depends on the model pa- 
rameters, and is known to be NP-hard when the interactions are sufficiently strong. Our 
first main result establishes that the computational transition for such models on (i-regular 
graphs is located precisely at the uniqueness threshold (see Defn. 1.6) for the corresponding 
model on the (i-regular tree. 

Theorem 1. For d > 3 and A > Xc{d) = (^Ld , unless NP = RP there exists no fpras for 
the partition function of the hard-core model with fugacity A on (i-regular graphs. 
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2 A. SLY AND N. SUN 

The transition point Xc{d) is the uniqueness threshold for the hard-core model on the 
d-regular tree: it marks the point above which distant boundary conditions have a non- 
vanishing influence on the spin at the root. In a seminal paper [18], Weitz used computational 
tree methods to provide a fptas for the partition function of the hard-core model on graphs 
of maximum degree d at any A < Xc{d). Together with Weitz's result, Thm. 1 completes the 
classification of the complexity of the hard-core model except at the threshold Ac. 

Previously it was shown that there is no fpras for the hard-core model at Xd > 10000 [13]. 
In the case of A = 1 this was improved to d > 25 [6, 7], using random regular bipartite 
graphs as basic gadgets in a hardness reduction. Mossel et al. [15] showed that local MCMC 
algorithms are exponentially slow for A > Xc{d), and conjectured that Ac is in fact the 
threshold for existence of an fpras. 

The first rigorous result establishing a computational transition at the uniqueness thresh- 
old appeared in [17], where hardness was shown for Xc{d) < A < Xc{d) + e{d) for some 
e(d) > 0. The proof relies on a detailed analysis of the hard-core model on random bipartite 
graphs, which are then used in a randomized reduction to max-cut. More precisely the re- 
sult of [17] gives hardness subject to a technical condition which was an artifact of a difficult 
second moment calculation from [15], and which could only be verified for A < Xc{d) + e{d). 
Hardness was subsequently shown by Galanis et al. [8] for all A > Xc{d) when d 7^ 4, 5 by 
verifying the technical condition of [17]. 

In this paper we follow a different approach which is more conceptual and completely 
circumvents second moment method calculations. Moreover the same method of proof gives 
the analogous result for anti-ferromagnetic Ising models with arbitrary external field: 

Theorem 2. For d > 3, B e R and (3 < /3c,af (-B, d) < 0, unless NP = RP there does not 
exist an fpras for the partition function of the anti-ferromagnetic Ising model with inverse 
temperature /3 and external field B on d-regular graphs. 

Here /3c,af(-B, c^) denotes the uniqueness threshold for the anti-ferromagnetic Ising model 
with external field B on the d-regular tree. Extending the methods of Weitz [18], Sinclair 
et al. [16] (see also [12]) gave a fptas for the anti-ferromagnetic Ising model on d-regular 
graphs at inverse temperature (3 > /3c,af (-B, d), so together with Thm. 2 this again establishes 
that the computational transition coincides with the tree uniqueness threshold. 

The hard-core and anti-ferromagnetic Ising models together encompass all (non-degenerate) 
homogeneous two-spin systems on d-regular graphs (see §2.2). Thus, the results of [18, 11, 16] 
combined with Thms. 1 and 2 give a full classification of the computational complexity of ap- 
proximating the partition function for (homogeneous) two-spin systems on d-regular graphs, 
except at the uniqueness thresholds Xc{d) and /3c,a.i{B,d). 

In fact, we will show inapproximability in non- uniqueness regimes in a strong sense: not 
only does there not exist an fpras, but for any fixed choice of model parameters and d 
there exists c > such that it is NP-hard even to approximate the partition function within 
a factor of e'^" on the class of rf-regular graphs. 

Independent results of Galanis— Stefankovic—Vigoda. In a simultaneous and indepen- 
dent work, Galanis, Stefankovic and Vigoda [9] established the result of Thm. 1, and Thm. 2 
in the case of zero external field {B = 0). Their methods differ from ours: they analyze the 
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second moment of the partition function on random bipartite d-regular graphs, and estab- 
hsh the condition necessary to apply the approach of [17]. Their proof analyzes a difficult 
optimization of a real function in several variables by relating the problem to certain tree 
recursions. 

1.1. Reduction to max-cut via bipartite graphs. Our proof is based on a detailed 
characterization (Thm. 5) of the local structure of anti-ferromagnetic two-spin systems on 
symmetric bipartite d-regular locally tree-like graphs. Specifically, we show that the joint 
distribution of all the spins in a large neighborhood of a typical vertex in the graph converges 
to a known (Gibbs) measure on the (i-regular tree. Under the additional assumption that 
that the graph is an edge expander, when the model has non-uniqueness on the (i-regular 
tree the spin distribution on the graph is divided into + and — phases where one or the 
other side of the graph has a linear number more vertices with + spin. 

Our main results Thms. 1 and 2 are then proved by a variation on the construction of [17], 
using the bipartite graphs in a randomized reduction approximate max-cut on 3-regular 
graphs, which is known to be NP-hard [1]. First, we use Thm. 5 to construct a symmetric 
bipartite (i-regular locally tree- like graph G of large constant size such that, conditioned on 
the phase of the global configuration, spins at distant vertices are asymptotically independent 
with known marginals depending only on the side of the graph (Propn. 4.2). 

Given a 3-regular graph H on which we wish to approximate max-cut, first we take a 
disjoint copy G^ of G for each vertex v & H . After removing 3fc edges from each Gj,, for 
each edge (u, f ) G H we add k edges joining each side of Gu to the corresponding side of G^, 
in such a way that the resulting graph H'~^ is (i-regular. 

The connections between gadgets do not substantially change the spin distributions inside 
them, and in particular the ± phases remain. The anti- ferromagnetic nature of the interation, 
however, results in neighboring copies of G in H^ preferring to be in opposing phases. Using 
the asymptotic conditional independence result Propn. 4.2 we can estimate the partition 
function for the model on H'^ restricted to configurations of given phase on each copy 
of G within a factor of e^'^' (Lem. 4.3). We find that the distribution is concentrated on 
configurations where the vector of phases gives a good cut of if, and the effect is strengthened 
as k is increased. Thus, for any e > 0, by taking k (hence G) to be sufficiently large a (1 + e)- 
approximation of max- CUT (if) can be determined from the partition function of the model 
on if '^, thereby completing the reduction. 

Our reduction depends crucially on the detailed picture of the spin distribution developed 
in Thm. 5 and Propn. 4.2. Using methods developed in [14], these results in turn are 
obtained as consequences of precise asymptotics for the partition function of two-spin models 
on bipartite (i-regular graphs: we show that the log-partition function, normalized by the 
number of vertices in the graph, has an asymptotic value, the "free energy density," which is 
easily computed from the (non-rigorous) "Bethe prediction" of statistical physics (see §2.1). 
This is a result of independent interest, since lower bounds for partition functions on graphs 
have proved to be in general challenging. Asymptotics for the partition function on general 
tree-like graphs were established for the ferromagnetic Ising model in [3, 5, 4], and for more 
general spin systems in uniqueness regimes in [4]. Our result for anti-ferromagnetic models 
is stated somewhat informally as follows; for the precise statement see Thm. 4. 
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Theorem 3. For any non-degenerate homogeneous two-spin model on bipartite (i-regular 
locally tree-like graphs, the log-partition function normalized by the number of vertices has 
an asymptotic value which coincides with the Bethe free energy prediction. 

In the remainder of this introductory section we formally introduce the models which we 
consider. We then define the notion of local (weak) convergence of graphs and give precise 
statements of our results on the partition function (Thm. 4) and local structure (Thm. 5) of 
these models on bipartite graphs. 

1.2. Definition of spin systems. Let G = {V, E) be a finite undirected graph, and ^ a 
finite alphabet of spins. A spin system or spin model on G is a probability measure on the 
space of (spin) configurations a_ E ^^ of form 

— {iJ)<^E i€V 

where ip is a. symmetric function J?r^ — )■ M>o, ip is a. positive function ^ — t- M>o, and Zci'ip) 
is the normalizing constant, called the partition function. The pair ip = {ip,ip) is called a 
specification for the spin system (1.1). 

In this paper we consider spin systems with an alphabet of size two; without loss ^ = 
{±1}. The Ising model on G at inverse temperature /3 and external field B is given by 



The hard-core (or independent set) model on G at activity or fugacity A is given by 

-^(-)=^ n i{-^-.^i}n^'' (1-3) 

where a = l{a = +1} = (1 + cr)/2. The edge interaction has no temperature parameter 
and includes a hard constraint. Our definition (1.3) is trivially equivalent to the standard 
definition of the hard-core model which has spin in place of —1, but we take ^ = {±1} 
throughout to unify the notation. 

1.3. Local convergence and the Bethe prediction. If G is any graph and v a vertex 
in G, write Bt{v) for the subgraph induced by the vertices of G at graph distance at most 
t from V, and dv = Bi{v)\{v} for the neighbors of v. We let T = {T,o) denote a general 
tree with root o, with T* = Bt{o) C T the subtree of depth t. We also fix d throughout and 
write T = (T, a) for the rooted d-regular tree. 

Definition 1.1. Let Gn = (Vn = [n], En) be a sequence of (random) finite undirected graphs, 
and let /„ denote a uniformly random vertex in Vn. The sequence Gn is said to converge 
locally to the li-regular tree T if for all t > 0, Bt{In) converges to T* in distribution with 
respect to the joint law P„ of {Gn, In)'- that is, lim„_^oo IPn(-Bt(-^n) — T*) = 1 (where = denotes 
graph isomorphism). 

We write E„ for expectation with respect to P„ and impose the following integrability 
condition on the degree of /„: 
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Definition 1.2. The sequence Gn is uniformly sparse if the random variables \dln\ are 
uniformly integrable, that is, if 

lim limsupE„[|a/„|l{|a4| > L}] = 0. 

We assume throughout that Gn (n > 1) is a uniformly sparse graph sequence converging 
locally to the c?-regular tree T; this setting is hereafter denoted Gn -^loc T. The free energy 
density for a specification ip on G„ is defined by 

= lim 0„ = lim -E„[logZ„], Z„ = Zg„{iP), (1.4) 

n— >oo n— >oo n — 

provided the limit exists. For ferromagnetic spin systems on a broad class of locally tree-like 
graphs, heuristic methods from statistical physics yield an explicit (conjectural) formula for 
the value of 0, the so-called "Bethe prediction" $ whose definition we recall in §2.1. For anti- 
ferromagnetic two-spin models, the Bethe prediction is well-defined only on graph sequences 
Gn which are nearly bipartite, in the following sense: let T+ denote the d-regular tree T 
with vertices colored -|-1 (black) or —1 (white) according to whether they are at even or odd 
distance from the root o; let T_ be T+ with the colors reversed. Let T be the random tree 
which equals T+ or T_ with equal probability; write P for the law of T and E for expectation 
with respect to P. 

Definition 1.3. For G„ —^loc T, we say the Gn are nearly bipartite, and write Gn -^loc T 
(equivalently Gn -^loc P), if there exists a (not necessarily proper) black-white coloring of 
Gn such that for all t > 0, Bt{In) — > T* in distribution. 

The precise statement of Thm. 3 is then as follows: 

Theorem 4. Let ip specify a non-degenerate homogeneous two-spin system. 

(a) If ip is ferromagnetic, then exists for any Gn -^hc T and equals ${t} as defined by 
(2.2) (and given more explicitly by (2.4)). 

(b) If ip is anti-ferromagnetic, then exists for any G„ -^loc T and equals ${t±} as defined 
in (2.2) (and given more explicitly by (2.3)). 

Remark 1.4. Hereafter we treat Gn -^loc T and Gn -^loc T in a unified manner when 
possible by writing Gn -^loc ^r fo^^ ^T the uniform measure on T, which always denotes 
either {T} or {T±}. We write E7- for expectation with respect to P7-. 

1.4. Local structure of measures. Under some additional assumptions on G„, Thm. 4, 
together with the arguments of [14], characterizes the asymptotic local structure of the spin 
systems z/„ = uq^- For Gn -^loc T, let r : V„ — )■ {±} denote the given black-white coloring of 
the vertices of Gn (hereafter writing ± as shorthand for ±1). We say that Gn is symmetric 
if it is isomorphism-invariant to reversing the black-white coloring. For a spin configuration 
a G Gn we define the phase of a to be 

Y{a) = sgnN^Tjai, where sgnx = l{x > 0} — l{x < 0}. 
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Let z/^ denote the measure z/„ conditioned on the configurations of ± phase: that is, 

" {ij)€E„ i&V„ 

where Z^ is the partition function restricted to the it configurations. We will characterize 
the local structure of the measures v^ on graph sequences satisfying an edge-expansion 
assumption, as follows: 

Definition 1.5. A graph G = {V,E) is a {S,'y,X)-edge expander if, for any set of vertices 
S C V with 6\V\ < \S\ < 'y\V\, there are at least X\S\ edges joining S to V\S. 

The measures z/^ will be related to Gibbs measures on the infinite tree. In particular, 
recall the definition of (Gibbs) uniqueness: 

Definition 1.6. For a rooted tree T, let ^t denote the set of Gibbs measures for the 
specification ip on T. The specification is said to have (Gibbs) uniqueness (on T) if \'^t\ = 1- 

Recalling Rmk. 1.4, let ^7- denote the space of mappings z/ : T h-)- z^(T), T E T (with 
^{T} ^^ ^{T±} in the obvious manner). When T = {T-t} we write i'± as shorthand for i'{T±). 

Definition 1.7. An element u G ^r is translation-invariant if for (T, o) E T and any vertex 
X E T, the law on spin configurations of (T, x) induced by z/(T, o) coincides with z/(T, x)} 

For a two-spin model, let z/"*" (resp. z/~) be the elements of ^7- defined by conditioning on 
all spins identically equal to 1 on the t-th level of black (resp. white) vertices and taking the 
weak limit as t — > 00; the z/^ are translation- invariant. The projections fi'^ = u^ = z/+(T+) 
and fi~ = u^ = z/~(T_|_), disregarding the black- white coloring on T+, are the extremal 
semi-translation-invariant Gibbs measures for the model on T, and by symmetry 

j2+ = uZ = Jy'{T-), ;u" = Z/+ = z/+(T_). 

The model has uniqueness if and only if fi^ = fi~. 

Definition 1.8. For Gn ~ Pn a random graph sequence and Un any law on spin configurations 
a^ of Gn, we say that P„ ® z/„ converges locally (weakly) to P7- (g) u (for z/ G ^r); and write 
Pn ® z/„ -^loc Pr ® z^, if it holds for allt > that {Bt{In),o_Q^^^)) converges in distribution 
to (T*, aj where T ~ P7- and a^ is the restriction to T* of a ~ i^{T). 

Remark 1.9. In [14, Defn. 2.3] three forms A,B,G of local convergence of measures are 
distinguished, with G ^ B ^ A. Our Defn. 1.8 corresponds to the weakest form A: however, 
as explained in the proof of [14, Thm. 2.4 (II)], if the {v{T))xeT are extremal Gibbs measures 
then A,B^G are easily seen to be equivalent, so convergence in the sense of Defn. 1.8 implies 
convergence in the a priori stronger sense of 

i|Pn[(5*(/n),^B.(/„)) = ■] -Pr[(T*,a^O = ■]||xv -^ 0. 
Theorem 5. For any anti-ferromagnetic two-spin system on Gn -^loc T, the following hold: 
(a) If the Gn are symmetric, then P„ (g) z/„ -^loc P ® [{^^ + '^ )/2]. 



If T — {T} this agrees with the usual definition of translation-invariance, whereas if T = {Tj-} then the 
projections ;^(T ) are serrei-translation- invariant. 
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(b) If for all (5 > the Gn are (5, 1/2, A5)-edge expanders for some A^ > 0, then 

P„, ® l^n ^loc P ® Jy^- (1-5) 

Further, with ( )^ denoting expectation with respect to the Gibbs measure /i, 

-y{^)'^Tiai -^ 2^^°)^*+ ~ (o"o)|i-] in probability. (1.6) 

Outline of the paper. In §2 we review the Bethe prediction in the (i-regular setting and 
prove Thm. 3 (in its form Thm. 4). In §3 we show how to deduce Thm. 5 from Thm. 4 by 
the methods of [14]. In §4 we prove the approximate conditional independence statement 
(Propn. 4.2) and demonstrate the randomized reduction to MAX- CUT to prove our main 
results Thms. 1 and 2. 

2. Partition function for two-spin models 

In this section we prove Thm. 4, establishing the free energy density (and verifying the 
Bethe prediction) for two-spin models on graph sequences Gn -^loc T. We refer to [2, 4] for 
more general background and references on the Bethe prediction, and in §2.1 describe only 
its specialization to the (i-regular setting. In §2.2 we show that for purposes of computing 
on (i-regular locally tree-like graph sequences, all non-degenerate two-spin systems reduce to 
Ising or hard-core. In §2.3 we compute the free energy density for these models by applying 
an interpolation scheme described in [4], thereby completing the proof of Thm. 4. 

2.1. The Bethe prediction. Recalling the notation of Rmk. 1.4, we now review the Bethe 
prediction for G„ -^loc Pr- Given T, let Te denote the set of trees T rooted not at a vertex 
but at an oriented edge x ^ y, obtained by distinguishing an oriented edge in T G T and 
forgetting the root. Elements of T, Te are regarded modulo isomorphism: thus if T = {T} 
then T; = {(T, o ^ j)}, and if T = {T±} then Te = {(T±, o -^ j)}. 

Let A denote the (|^| — l)-dimensional simplex of probability measures on <^. A message 
is a mapping /i : 7^ — )■ A; we write Ti = 'H{T) for the space of messages on %. For T E T, 
a; — )■ y in T, and h ETi, write h^^y for the image of (T,x ^ y) E % under h, and define 

where 

«l>^^(/i) ^log|^7/;(or„) II (j2ijiao,aj)hj^,iaA\, 

^ (Jo j€do ^ aj ^ ^ 

jedo ^ o-o,a-j 

The Bethe free energy functional on 'H(T) is defined by $r(^) = ^Ti^rih)]. 
The Bethe or belief propagation (BP) recursion is the map 

BP = BPr : H{T) ^ H{T), {BPh),^y{a) = F[(/i,^J,ga.\,] 
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for F : A''"^ — > A defined by 

[F(^)](cT) = V^(a)n|E^(^'^^)^^(^^)|' h^{hi,...,h,.,)eA'-' (2.1) 

(where = denotes equivalence up to a positive normalizing factor). 

Definition 2.1. For any homogeneous spin system on G„ -^loc Pr? the Bethe prediction is 
that the free energy density (p of (1.4) exists and equals 

$ = $r = sup ^r{h) (2.2) 

hen* 

with n^ = n^iT) C n{r) the set of all fixed points of BPr- 

For /i G A write F(/i) = F{h, . . . ,h): then "H^dT}) corresponds simply to the fixed points 
of F in simplex. For h G 'H{{T±}) we write h±. = h(T±, o -^ j) G A: then any h G ^^^,{{Tj-}) 
must satisfy h± = F(/i=p), so ^{^{{T^}) corresponds to the fixed points of the double recursion 
F(^)^FoF. 

In verifying the Bethe prediction we will identify the fixed points attaining the supremum 
in (2.2). In the anti-ferromagnetic case, with h~^ (resp. h~) denoting the elements h G 
'Hi,{{T±}) maximizing /i+(+) (resp. /i_(+)), we will see that 

${T±} = ${T±}(/^+) = ${T±}(/^~). (2.3) 

Explicitly, h'^ = hZ (resp. /i^ = hl^) will be the fixed points of F'^-* giving maximal (resp. 
minimal) probability to spin +. The ferromagnetic case reduces to the Ising model: here, 
with h^ denoting the elements of 'H^,{{T}) maximizing ho^j{±) on T, we will see that 

${T} = $m(/^^^"''). (2.4) 

The remainder of this section is devoted to the proof of Thm. 4. 

2.2. Reduction to Ising and hard-core on rf-regular graphs. We first show that for the 
computation of the free energy density, all (non-degenerate) homogeneous two-spin models 
on graph sequences G„ -^loc T reduce to either the Ising or hard-core model. Indeed, let 
ip = (ipyi^) be a specification for a two-spin system with alphabet ^ = {±}. If we define 
^' by i/''((7,a') = ^lj{a,a')i>{ay/'^i>{a'y/'^, and i>'{a) = 1, then 

^logZcif) -UogZcif) = 0(E„\\dIn\l{\dIn\ ^ d}]) , 

which for G„ -^loc T tends to zero as n — ;■ oo by uniform sparsity. Therefore we assume 
without loss ip = 1, and consider the possibilities for tp: 

(1) If z^ > 0, then V^(a, a') = e^^e^''''' e^^/'^e^"'/'^ for f3, B, Bq defined by 

SO (pn — {d/2)BQ is asymptotically equal to the free energy density for the Ising model 
on Gn with parameters {/3,B). 



^0 "^' "'^ +B + -En 

n n 
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(2) IfV^(+,-) = V^(-,+) > OandV^(-,-) > ip{+,+) = 0, then, recalling ct = l{a = +}, 
we have ^p{a, a') = e^n{aa' ^ IjX^/'^X^'/'^ for Bq, A defined by 

W[ , )-e , ^(_^_) -^ • 

Therefore 0„ — {d/2)BQ is asymptotically equal to the free energy density for the 

independent set model on G„ at fugacity A. 
The remaining two-spin models are degenerate, with free energy density which is easy to 
calculate: 

(3) Suppose tp{+, — ) = ^{—, +) = 0, so that ■?/'((T, a') may be written as l{a = Q-'JQ^OQBa/d^Ba'/d^ 

Then 

■HG„) -. 

5^1og(l + e-2^l^^l) 
• i=i -' 

where the sum is taken over the connected components Ci, . . . ,Ck{G„) of Gn- We 
claim 0„ — )■ = {d/2)BQ + B: we have liminf„^oo(0n — 0) > (using uniform 
sparsity), and 

hm sup {(pn — 0j < hm sup log 2 , 

ra— >oo n— >oo "^ 

SO it suffices to show E„[A;(G„)]/n — )• 0. Indeed, if this fails then there exists e > such 
that for infinitely many n, the event {k{Gn) > e^} occurs with P„-probability at least 
e. On this event. On has at least en/2 components of size < 2/e, so for t > logji.(2/e), 
limsup„_,^P„(5i(/„) ^ T*) > eV2 > 0, in contradiction of G„ ^loc T. 

(4) Suppose instead ip{+, +) = '*/'(—, — ) = while ■?/'(+, — ) = '?/'(—, +) > 0. If the G„ are 
not exactly bipartite then 0„ = — oo. If they are exactly bipartite then 

0„ = log^(+,-) + log2 , 

n n 

and by the observation of (3) this converges to = {d/2) log ■?/'(+, — )• 

2.3. Bethe interpolation. We now evaluate the hard-core and Ising free energy densities 
by interpolating in the model parameters. Write ^ = log-?/', ^ = log'0', and for the hard-core 
model take B = log A. Let ( )^'^ denote expectation with respect to v!^'^ = Vq , and define 

ar(/3,5) = 5b0„(/3,5) = M{dBi{cri^)f^% 
(with a^{P,B) = for hard-core). We also define analogous quantities on the limiting tree 



for 

T ~ Pr: ioT hen let 



a™(/3, B, h) = af^(/3, 5, h) = Er[ldBa^o)t'% 
a\P, B, h) ^ a^(/3, B, h) ^l^r[j2 1^/3^(^- ^^•)1^'1 (2.5) 



jGdo 
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where {-Jh denotes expectation with respect to the measure i/^i on spin configurations on T^ 
defined by 

jedo 

The following lemma, describing our interpolation scheme, may be verified directly or ob- 
tained as a consequence of [4, Propn. 2.4]. We always interpolate in one parameter at a time, 
keeping the other fixed and suppressing it from the notation. 

Lemma 2.2. // for B G [-Bq? -Bi] we have h = h{B) G "H^ which is continuous and of 
hounded total variation in B, then 

$r(5i) - $r(fio) = / af^(5, h) dB. 

J Bo 

The same result holds for B, a™ replaced with 13, a^. 

We now make explicit the connection between BP fixed points and Gibbs measures; for a 
discussion in a more general setting and further references see [4, Rmk. 2.6]. Recall that for 
T E T, ^T denotes the set of Gibbs measures for the specification ip on T, and Wj- denotes 
the space of mappings T h-> z/(T) G ^t with T eT. An element /i = v{T) G ^t is a Markov 
chain or splitting Gibbs measure (see [19]) if there exists a collection h^^ = {h'^_^y) of elements 
of A indexed by the oriented edges of T such that for any finite connected induced subgraph 
U={Vu,Eu)o{T, 

where p{j) denotes the unique neighbor of j inside U for j belonging to the external bound- 
ary dU of U . Extremal Gibbs measures are Markov chains but the converse is false. We say 
that an element v G ^7- is Markovian if z^(T) is a Markov chain for each T E T: the associ- 
ated collection h" = (h''^^^)T£T is called an entrance law: it satisfies consistency conditions 
imposed by (2.6) (which closely resemble the BP equation), and the correspondence between 
Markovian z/ G ^7- and entrance laws h'^ is bijective. If u is also translation-invariant (in 
the sense of Defn. 1.7), then each hx-^y depends only on the isomorphism class of (T, x ^ y) 
in Te, so h E 'H(J~), and in fact by the consistency conditions h G 'Hi.iT). Thus there is a 
bijection between BP fixed points h G 'Hi.iT) and translation- invariant Markovian z/'* G ^r- 
In particular, for two-spin models, the z/^ G ^r of §1-^ and the h"^ G "H^ of §2.1 are related 
by this correspondence and so may be regarded as essentially equivalent. 

The main implication of Lem. 2.2 is the following (which may also be obtained as a special 
case of [4, Thm. 1.13]): if for B G [-Bq, -Bi] we have h = h{B) G Ti^ which is continuous and 
of bounded total variation in B, then 

limsupa™(5) < a^^{B,h), (2.7) 
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implies limsup„^oj3[0„(i?i) — ^^(-Bo)] < *^r(-Si) ~ ^t{Bo). (The statement also holds with 
B, a™ replaced by /3, a^.) From the above discussion we can re-express 

jedo 

so showing (2.7) amounts to proving a relation between the expectation of a local observable, 
in this case 9b^((Jj), in the finite graph to the expectation of the analogous observable 
ds^io'o) under translation- invariant Markov chains on the limiting tree. In the remainder of 
the section we carry out this scheme. Note that for the Ising model 

aniB) = E4{ajX], <(/?) = l^n[j2 (^^"^^■)n" ■ 

For the hard-core model a^{B) = E„[((T/^)^]. 

2.3.1. Interpolation for hard-core. 

Lemma 2.3. For the hard-core model at fugacity X, the supremum of {cjo + d~^ ^^jedo^j)!^ 
over fi & '^j is achieved precisely by the measures fi^. 

Consequently, the supremum of E[{d'o)u] over translation-invariant v G ^r ^s achieved pre- 
cisely by the p^ . 

Proof. By extremal decomposition, assume without loss that /i is itself extremal, with h^ = 
{h'^-^y) as defined above. For j G do write qj = h^^^{—): then 

/ l>r- \ _ ^ rijgoo gj + d'^ Ejgao(l - (Ji) 

dl + \ rijeao (Ij ' 
For fixed ^ . q^ this is (strictly) maximized by taking all q^ = g, so the above is 

maxg-<,<g+ [g~i + Ag^-^] ' 

where q~ and g"^ are the minimal and maximal values for g, corresponding to /i~ and yU"*" 
respectively. Since g~^-|-Ag'^~^ is convex, the maximum can only be attained at the endpoints, 
and in fact it is attained at both endpoints with value (l/g~ + 1/g^ — 1)^^- D 

Lemma 2.4. For the hard-core model, "H^dT}) consists of a single message h* = h*{X). 
'H^{{T±}) consists of the messages h*,h^ which coincide for X < Xc and are distinct for 
X > Xc. The messages are continuous in X, smooth except possibly at X = Xc- 

Proof. For the hard-core model, the function F(/i) {h G A) of (2.1) is expressed in terms of 
g = h{—) as F(g) = (l + Ag*^"^)"^. As g increases from to 1, F decreases from 1 to (l + A)"^, 
so F has a unique fixed point g^, which is smoothly decreasing in A. We compute 

F' = -^^F[l - F], F" = -^F(l - F)[{d - 1)2F - d], 
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SO G = F*- ^ has second derivative 

G" = (F' o F)F" + (F" o F)(F')^ = ^^^^^(1 - F)(l - G)G ■ Q, 

Q = -2{d - 1)2(1 - F)G +{d- 2)[d{l - F) + F]. 
Setting this to zero gives 



2{d-iy V i-F. 

The left-hand side is decreasing in F while the right-hand side is increasing, so G has at most 
one inflection point, hence at most three fixed points. If G has a fixed point which is not 
equal to g^,, then necessarily it has exactly three fixed points F(go) < q* < Qo, so G'(g^) > 1. 
But 

G'{q^) = F{q^y = {d-ir{l-q^y 

is smoothly increasing in A with G'(g^) = 1 precisely at A = Xc{d), so we see G has a unique 
fixed point /i^ when A < Ac, and when A > Ac it has three fixed points F(go) < q* < qo which 
are smooth on the open interval (Ac, oo). 

It remains to verify that go — > Q'* as A J, Ac- Suppose otherwise, so that 

limsup F(go) + 2e < g^, < liminf go — 2e 

for some e > 0. It is possible to take a sequence A J, -^c along which the inflection point of 
G always lies on the same side of g^: assume it is < g^, (the argument for the > g^ case is 
symmetric), so that G' is decreasing on g > g^. By the mean value theorem applied to the 
interval [g^ + e, g^, + 2e], 

G'(g. + e)> G(^- + 2e)-G(g. + e) ^ q. + 2e - [q. + G\q.)e] ^ ^ _ ^,^^^^_ 

e e 

Therefore 2 — G'(g^) < G'(g) < G'(g^) for all g G [g*,g* + e], and consequently F'(g) = 1 for 
all g G [g*, q^, + e] at Ac- This gives the desired contradiction and the lemma follows. D 

Proposition 2.5. For the hard-core model, 

(a) The definitions (2.2) and (2. -3) of ^ = ${t±} coincide. 

If Gn -^loc T then = $ for A < Ac and lim sup^ 0n < $ for A > Ac- 
(h) IfGn ^loc T then <p = ^ for all A > 0. 

Proof, (a) Since any subsequential local weak limit of the measures Vn must be translation- 
invariant, the second part of Lena. 2.3 implies 

limsup<^(E) = limsupE„[(a,jf] < E[(a„),±] = a'^{B,h^) 

n— >oo n— >oo 

(the inequality can alternatively be obtained by expressing a^{B) as o(l) + E„[(o-/^ + 
d^^J2jediu^i^n]/'^ sind directly applying the first part of Lem. 2.3). Also, a^^{B,h*) < 
a^^{B,h~^) = a^^{B,h~), with equality for A < Ac and with strict inequality for A > Ac. It 
then follows from Lem. 2.2 (using Lem. 2.4) that for Ac < A, 

limsup[0„(A) - 0„(Ac)] < <I>(A, h^) - $(Ac, h^) > $(A, h*) - $(Ac, h*). 
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It was shown in [4, Thm. 1.11] that = $(/i*) = $(/i^) for A < Ac so the claim follows 
(again making use of Lem. 2.4). 

(b) It suffices to show that for Gn -^loc T, 



Indeed, 



lim liminf(0„ - $) > 0. (2.8) 

A— ^oo n^oo 



<|.-(A) = ^ log[A(g+)'^ + 1] + ^ log[A(g-)'^ + 1] 



$<^(A) = ^log(l-(l-g-)(l-g+)), 



and limA^oo q~^ = hmA^oo(l - g ) = 1, so limA^oo $(A) - [log(A + l)]/2 = 0. But 



liminf (/)„ > liminf — E 



n 
n— >oo n— >-oo TL 






log(l + A)- 



n 



for an the independence number of Gn- But an is at least the number of black vertices with 
no black neighbors, so Gn -^loc T implies limsup„E„[a„]/n < 1/2. This proves (2.8) from 
which the result follows. D 

2.3.2. Interpolation for Ising. 

Lemma 2.6. For the Ising model with parameters /3 < and B E M., the supremum of 
{ao + d^^ Xljeao ^i)^ ^^^''^ fi E ^j is achieved precisely by the /x^. 

Consequently, the supremum ofE[{ao)u] over translation-invariant v G ^j is achieved pre- 
cisely by the z/^. 

Proof. We argue as in the proof of Lem. 2.3: assume fi is extremal, and write {hj, 1 — hj) = 
(/^;^.(+),/^^_,,(-)): then 

where R^ = [e-^hj + e^(l - hj)]/[e^hj + e'^{l - h^)], and 

^^e%-e-P{l-h,] 



B. 



e^hj + e-/5(l - hj) 
e-^h, - e^{l - h,] 



2Rj 


— 


(e 


2/^ + e 


-2.) 


{e'^ 


e~ 
+ 


-2/3 

e' 


_e2/3 

-2/3) _ 


5 



' ~ e-Phj + e/5(l - hj) e-2/3 - e^/^ 

If {hk)kj^j are fixed, the expression is maximized by taking hj as large or as small as possible. 
Since /3 < 0, for any fixed YlkedoR-k, both J2jedo^j ^^^ J^j^do^j ^^^ maximized by taking 
all the Rj equal, i.e. with hj = h. The overall maximum is then attained for h equal to h~^ 
or h~ , the minimal and maximal values for q corresponding to the /i^. In fact it is attained 
at both endpoints with value 

2{h- + h+ -1) 
1 + (e-2/3 - l){h- + h+) - 2(6-2/3 - l)h-h+' 



14 A. SLY AND N. SUN 

which concludes the proof. D 

Lemma 2.7. For the Ising model with parameters /3 < and B = 0, the supremum of 
~{J2jedo^oO'j)fi over fi E ^j is achieved precisely by the /x^. 

Proof. As in the proof of Lem. 2.6, assume /i is extremal and write {hx^y, 1 — hx^y) = 
{K^y{+),K^y{-)). Then 

_, . _ e-^ - (e^ + e-^)[K^,hj_, + (1 - /i,^j)(l - /i,^,)] 

^""^ '^^ e-/3 - (e-/^ - e/^)[/i._,/i,^o + (1 - /io^,)(l - h,^o)] ' ^ ' ' 

The partial derivative with respect to hj-,o has the same sign as 1/2 — ho^j, so —{cToO-j)^ is 
maximized with {ho-^j, ho^j) equal to (/i^, /i^) and {h~^, h~), corresponding to the measures 

Lemma 2.8. For the anti-ferromagnetic Ising model, "Hf'^dT}) consists of a single message 
h* = h*{f3,B). 'H^{{T±}) consists of the messages h*,h'^ which coincide for (3c,a.fiB,d) < 
/3 < and are distinct for [3 < f3c,i^f{B,d). The messages are continuous in /3,B, smooth 
except possibly where (3 = Pcm{.B, d). 

Proof. For the Ising model, the function F(/i) {h G A) of (2.1) is expressed in terms of 
t = \og[h{+)/h{-)] as 

F(t) = 2i? + (d-l)log|^^|, 9^e-'^ (2.10) 

F — 2B is an odd function of t G M, identically zero when /3 = and strictly monotone 
otherwise, going from {d — 1) log 6* to —{d — 1) log 6' as t increases from — oo to cx3. Suppose 
from now on that /3 < 0. Then F has a unique fixed point t* of the same sign as B, smoothly 
increasing in B, and smooth in (3 with absolute value increasing as (3 becomes more negative. 
For G = F^^^ we compute 

^'~ (eF + ^)2(^eF + l)2(e* + ^)2(^e* + l)2 ' 

A± = -e{e^' - 1) ± (rf - l){e^ - l)e*, 

A = -{e^ + l){e^'-l). 

If A+ = then clearly G can have at most one inflection point, so suppose A+ > 0: then 
setting G" to zero results in e^ = r± where 

r±^(-A±v^)/(2A+), 

D = A''- AA+A^ = {6^ - 1)2[1 + 4((rf - 1)^ - 1/2)6'* + e^*] > 0. 

If A/A^ > then at most one of the r±. can be positive. If A > then t < which implies 
A^ > 0, so it remains only to consider the case A < < A+: in this case, A < implies 
t > and so A_ < 0, therefore 

\D\ -A^ = -4:A+A_ > 

which implies r_ < < r_|_. Thus G has at most one inflection point for any /3 < 0. 
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By implicit differentiation we find 



[dom]\t=t. 



1 - F'(t 
(rf - l)e* [die^ - 1) + 2] (e^* + 1) + 4^e* 



<0, 

t=t^, 



liminf<(/3) = -liminfE„ ^^ (cr/„cri)( 



(e* + ^)(^e* + 1) [d{e^ - 1) + 2]e* + ^(e^* + 1) 

so G'(t^) = F'(t^)^ increases smoothly as f3 decreases. The result then follows by repeating 
the argument of Lem. 2.4. D 

Proposition 2.9. For the anti-ferromagnetic Ising model, 

(a) The definitions (2.2) and (2.3) of ^ = ${t±} coincide. 

If Gn -^loc T then (jy = ^ for ^ > /3e,af(5) and limsup„ 0n < $ /or /3 < /3c,af(5)- 
^5; // G„ ^loc T t/ien = $ /or a// /3, E. 

Proof, (a) First fix i? = 0: since any subsequential local weak limit of the measures Un must 
be translation-invariant, Lem. 2.7 gives 

Also, a'^(/3, /i*) > a'^{f3, h^) with equality for [3 > /3c,af and with strict inequality for (3 < /3c,af- 
It then follows from Lem. 2.2 (together with Lem. 2.8 and the previous result for /3 > 0) 
that for i? = and /3 < we have 

limsup(/)„ <<l>(/i^) > ^(/i*). (2.11) 

Using Lem. 2.6 to interpolate in B (as in the proof of Propn. 2.5) then gives (2.11) for all 

/3<0,5gM. 

(b) Consider the limits /3 — )■ — 00 and 5 — )■ 00: 

lim [$(/3, 0) + pd/2\ = 0, lim [$(/3, B) - B - pd/2] = 0. 

/3— >■— 00 _B— >oo 

If G„ ^locT, then 

liminf 0, > liminf-/3 ^"^'^"'^ = -/3rf/2, 

liminf 0„ >B + liminf /3^^^^^ = fi + /3c//2, 

SO in fact = $ for all /3, S G M. D 

For completeness we review what is known for the ferromagnetic Ising model: 

Proposition 2.10. For the ferromagnetic Ising model on Gn -^loc T, exists and equals 
${T} as defined by (2.2) (and given more explicitly by (2.4)^. 
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Proof. In this setting "Hj, corresponds simply to the fixed points of a single iteration of the 
map F of (2.10), which is analyzed for example in [4, Lem. 4.6]. F always has between one 
and three fixed points, and we write t~^ and t~ for the maximal and minimal fixed points 
respectively, corresponding to h~^ and h~ as in (2.4). By symmetry we may always suppose 
B>0. 

At -B = 0, a fixed point is always given by t° = 0, unique provided F'(0) < 1. F'(0) 
increases monotonically in /3 and reaches 1 at (ic,i{B = 0,d), and for /3 > /3c,f(0, d) there are 
three distinct fixed points t~ < = t° < t~^, with t"*" = —t~ J, = t° as /3 J, /3c,f(0, d). Since 
adding B simply shifts the map F of (2.10) by the constant 2B, it is easy to deduce the 
behavior for general B > 0: if f'{t)\t=o < 1 then F has a unique fixed point t~^ = t~ which 
is zero when B = Q and increases smoothly in 5. If F'(t)|t=o > I5 then at S = the map 
F has three fixed points t~ < t° = < t"*". As i? increases, t^ increase smoothly while t° 
decreases smoothly. The fixed points t^ and t° merge at the threshold B = Bc^i{(3,d), and 
for B above this threshold we again have t~ = t"*". 

It follows from [3, Thm. 2.4] (see also [4, Thm. 1.8]) that for Gn -^loc T, (p exists and 
equals $ as defined by (2.4). Therefore it remains to verify that 

^{h+) > ^{h-) V ^{h°) for all B > 0. (2.12) 

At B = 0, it follows from the above that $(/i°) and $(/i+) = $(/i~) are continuous in /3. 
Writing h = h{+) and recalling (2.9), we compute 

OhK^oOj) [g2/3_2(e2/3_ !)/,(! _/,)]2' 

so {(Joaj)h+ > {<Jo<7j)h°, hence a^^P.k^) > a''{l3,h°), and then Lem. 2.2 gives $(/i=^) > $(/i°) 
at -B = 0. Next, for all B > 0, clearly a^{/3, /i+) > a'^{(3, h~)\/a'^{/3, h°), so another application 
of Lem. 2.2 gives (2.12) from which the proposition follows. D 

Proof of Thms. 3 and 4- Follows by combining the reduction of §2.2 with the results of 
Propns. 2.5, 2.9, and 2.10. D 

3. Local structure of measures 

In this section we show how Thm. 4 can be used to deduce Thm. 5 by straightforward 
modifications of the arguments of [14]. 

Proof of Thm. 5 (a). Observe that 9|0„ = n~^Er,[{S^)^ - {{S)^)^] where S is ^-^y^ ai for 
Ising and Yli^v ^i ^^"^ hard-core, so the 0„ are convex and hence so is the limit 0. Convex 
functions are absolutely continuous, so it holds for a.e. B that 0„, are differentiable in 
B with dB(t)n -^ dB4> = 9b$ (by Thm. 4). It follows from [4, Propn. 2.4] that ds^ = 
a^^{B, h^) = E[((To)i,+]. But for any subsequential local weak limit u of the z/„, we also have 

dB<t>n = ar(S) = ^nlicTlX] ^ E[(a„)J. 

Therefore E[(cro)z/] = E[(cro),y+], and it follows from Lem. 2.3 and Lem. 2.6 that t^ is a convex 
combination of the z/^. Since the Gn are symmetric, we must have z/ = (i/+ + z^~)/2. D 

We now analyze the conditional measures z/^, beginning with an easy observation: 
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Lemma 3.1. For anti-ferromagnetic two- spin models on Gn -^loc T, 



lim E„ 

n—^ca 



; { y^ T-jCyj = 



0. 



ieV„ 

Proof. For the Ising model see [14, Lein. 4.1]. For the hard-core model, let A^ denote the 
set of vertices i E Vn with B2{i) isomorphic to T^, the depth- two subtree of T+; then An is 
necessarily an independent set of black vertices. The probability that YlieA '^i^i ~ YlieA ^« 
takes value j, conditioned on all the spins {(fi)i^An^ is F{X = j) where X is a binomial 
random variable on N = \{i & An : (Tqi = 0}| number of trials with success probability 
A/(l + A), li N > en then P(X = j) = 0{l/y/en) uniformly in j (e.g. by the Berry-Esseen 
theorem). If A^ < en then J2iedA ^i — (l^nl ~ ^n)/d, so 

1^ 2^ _ 1 >v- / \An\/n-e \Vr\{AnUdA^)\ 1 ^^ 

n ^-^ n ^ — ^ n ^-^ a n n ^-^ 

iGVn i&Vn i&Vn i£Vn 

As n — 7- oo the right-hand side tends in probability to [— l/2-|-e((i+l)]/(i, which is negative for 
small e. Combining the above observations concludes the proof for the hard-core model. D 

In view of Lem. 3.1 we may without loss restrict attention to the measures z/^. Define the 
local functions (cf. [14, eq. (3.9)]) 

Fl indicates the vertices of G„ which are locally not in the + phase. 

Proof of Thm. 5 (b). We outline the steps of the proof of (1.5) following [14], describing 
modifications where needed. 

• Let z/* denote any subsequential local weak limit of the u^. Then u* G ^7- (see [14, 
Lem. 3.4]). By Lem. 3.1, z/+ has free energy density converging to 0, so the proof of 
Thm. 5 (a) implies that u* = {1 — q)v^ + qv~ for some q € [0, 1]. 

• By local weak convergence, lim„_^ooE„[(Fj^J„] = E[(F^)j^*]; further, if J„ denotes a 
uniformly random neighbor of J„, then 

lim M{nFl ^ Fl})n] = E[(1{F* ^ F]})A, j e do 

n— )-oo 

(cf. [14, Lem. 3.7]). 

• For the hard-core or anti-ferromagnetic Ising model in non-uniqueness regimes, there 
exists 6 > such that 

lim(F„*),+ = 0=l-lim(F„*),-, 

I— >oo t— >oo 

hm (F„* ^ i^*),+ = = hm (F,* ^ Fj),- 

(cf. [14, Lem. 3.8]). It follows that for sufficiently large t 

lim E4{Fl)n] >q-e, lim E„[(1{F„* ^ Fn)„] < e. 
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The argument of [14, Propn. 3.9] (using the edge-expansion hypothesis) now gives a contra- 
diction unless g = estabhshing (1.5). The proof of (1.6) then follows from applying the 
proof of [14, Thm. 2.5] to the bipartite case. D 

4. Computational hardness 

In this section we construct the bipartite expander gadgets to be used in the reduction 
to MAX-CUT (Lem. 4.1) and refine Thm. 5 to an approximate conditional independence 
statement for the gadgets (Propn. 4.2). We conclude with the proof of our main results 
Thms. 1 and 2. 

For any fixed positive integer k, G2n ^i^^ ^e a bipartite graph on 2n vertices with n even, 
defined as follows: 

• Let Hn be a graph on n vertices of maximum degree d, generated by the configuration 
model as follows: take a uniformly random matching m of [dn], and put an edge {i,j) 
in Hn for every edge {i',j') G m with i' G i+nZ, j' G j+nZ (self-loops and multi-edges 
allowed) . 

• Take G2n to be the bipartite double cover of if„: the two parts of G2n are («+)?=! and 
{i-)2=i, and we put two edges (z+,j_) and (j+,i_) in G2n for every edge (z, j) G Hn 
(multi-edges allowed). 

• Choose k vertices (i^)^=i uniformly at random from Hn, and for each i choose j^ G di^ 
uniformly at random. G'2„ is the simple bipartite graph formed by deleting the edges 
{i^j-,jt) from G2n and merging any remaining multi-edges in the graph into single 
edges. Write W^ = (4, jiJti and W = W+UW-. 

The graphs (j2n are d-regular with probability bounded away from zero as ra — > cxd (see e.g. 
[10, Ch. 9]). 

Lemma 4.1. Let k be fixed. For all 6 > there exists A^ > such that the G\n O'l"^ 
{6, 1/2, Xs)-edge expanders with high probability as n ^ oo. 

Proof. By stochastic domination we may assume d = 3. For S C i/„ with l^l = m, the 
probability that there are exactly j edges in Hn between S and its complement is 

p _, (3;)(3(np))^-,M3(^_^.)M3(„-^-,) 

where Ij^^ is the indicator that m — j is even, and M^ = (£ — 1)!! = 7r"-^/^r[(£ + l)/2]2^/^ is 
the number of matchings on [£] for £ even. By Stirling's approximation, if 5 < m/n <1 — 5 
and j = 'jn, then 

Pj,m = Ij,mexp I - n -H{m/n) - 7 log 7 + 05(7) + 05{n)j 

(where H{p) denotes the binary entropy function —plogp — {1 — p) log(l — p))- There are 
< e^^i'^/^) subsets of Hn of size m so there exists 75 > such that with probability at least 
^g-n_H'(<5)/4^ all subsets of Hn of size between 5n and (1 — 5)n have expansion at least 7^. 

We now show expansion for Gg^: since k does not change with n and the number of edges 
leaving any set of vertices decreases by at most a factor of 3 when multi-edges are merged 
into single edges, it suffices to show expansion for G2n- Let S± be subsets of the ± sides of 
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G2n such that S = S+U S_ has size < n. If the projection nS of 5* in iJ„ has size < (1 — 6)n, 
then S has expansion at least 7^/2. Suppose {ttSI > (1 — S)n: without loss l^+l > l^-l, so 
|7rS'+\7rS'_| > (1/2 — 6)n. If there are fewer than 7|S'| edges leaving S, then there must be at 
least 3(1/2 — 5)n — 7n edges between 7rS'+\7rS'_ and its complement in Hn- A similar analysis 
as above shows that for sufficiently small S there exists 75 > such that the probabihty G2n 
has such a set S" is < e""^'"^^^/*^, and this concludes the proof. D 

Recall that we use W"^ to denote the endpoints on the ± sides of the 2k edges deleted from 
G2n in the formation of G2n- Recall also the definitions of /i^ G ^j, and write h"^ = h'^^j G A. 
For h,h' e A define h ®^ h' e A.3K-2 by 

{h^^h'){a,a') 

for z{h ®^, h') the normalizing constant. 
Proposition 4.2. The conditional measure u ,. {a 



h{a)i;{a,a')h{a') 



z{h®^ h') 



(4.1) 



'-'2n 



converges to the product measure 



Q 



w 



a 



n ^^(^-) n ^^(^'^ 



vew+ 



w&W' 



Proof. Let Bf denote the union of the balls Bt{w) C Ggn over w G {ij.}^^^; assume that Bt 
is a disjoint union of graphs isomorphic to T* with internal boundary St = Bt\Bt-i, which 
is the case with high probability. For rj G J^'^* let 



t,e,r]V) 



c 



so that 



^G2„[(0-ii,0-,-^ 



= 


^GL 


{cr^i 


= ■ 


^S,= 


V) 


^ 


^GL 


(^ii 


= ■ 


^S, = 


v) 


) 


= ■ 


^sr- 


= 1 


= ^t,e,ri 


®, 



V' M/,r;- 



By Thm. 5, the conditional measures I'G^S—Btii'^ ) ~ ') converge to /i"*". But by (1.6), Y{q) 
agrees with Yt{a) = sgn^^gyv^^ tjctj with high probabihty, so that convergence also holds if 
we replace z/+^^ by VqIS') - '^Gini' I Yt{a) = ±). In particular, 

1 
\Stl,ac 



lim lim E 



2n 



Jtp 



M,f,o-o /t 



h^ ®^ h' 



'St G2„ 

On the other hand, it is easily seen that {h^^h'){l, 0) is maximized by taking h{l) and h'{0) 
as large as possible. But in the limit t ^ 00 the values ^fe „i^) , Cfe nW (with 77 arbitrary) 
are sandwiched between h'^(l), so it must be that 







lim lim E 

t— >oo n— >-oo 



2n 



(1I&-5, -^^IItv + IKm,.,^-^ IItv)^ 



'g 

'-'2n 



(4.2) 
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We now claim that (4.2) continues to hold after removal of the edges (ij_, j£;). Indeed, 



5Lfe.=!z) z^MzM i:,_^zti(i)zm ^^^^ 



(j2 



^P.{^s,=ri) Z^MZtM j:,'Z:M)ZUrr 



where 



Zoutiv) = ZG2n\Bt-i[{o:G^„\Bt-i ■ ^tU) = ± and g_s^ = r^}], 

Now note that for k bounded and t large we have ^m(??) ^ Z^j^{rj) uniformly over rj: for Ising 
interactions at non-zero temperature this is obvious, while for the hard-core model 

Ziniv) 



7k M-UV^- &.(1)Q,.(1)][1 - ^^Avi^KJ^)] 

which for t large is x 1 uniformly over r]. Since the ^^ ^ and Ctiv are ?7-measurable, it follows 
from (4.3) that (4.2) continues to hold with vp^ in place of z/^* . Since the spins (cr^)t„gvK 
are independent under z/^^ (■ | a^J, this further implies 



= lim lim E2„ 



Finally, by a similar argument as before lim„_>oo ^g^ (^(^l) = ^t(oi)) = 1, so (4.4) holds with 
i/i in place of v^l which gives the result. D 



We now demonstrate how to use Propn. 4.2 to establish a randomized reduction from 
approximating the partition function to the problem of approximate max-cut on 3-regular 
graphs, which is NP-hard [1]. 

Let if be a 3-regular graph on m vertices and construct the bipartite graph G = G^ by 
the procedure described above. By Lem. 3.1 and Propn. 4.2, for any e > there exists n(e) 
large enough such that the following hold with positive probability: 

(l) Ggn '^as formed by removing 3fc distinct edges from a d-regular graph 6*2™; 
(II) z/G3^(>^(^) = +)<(l + e)/2;and 

(ill) i^^iJyq_^r)/Q%{^w) ^ [1 - e> 1 + e] for all £vy 

Consequently, for given e we may find G^ satisfying properties (l)-(lll) within finite time 
by deterministic search. We then construct from H and G a new graph H^ as follows: 

• For each vertex x G ii let Gx be a copy of G, and denote by W^ the vertices of Gx 
corresponding to W^ in G. Let H^ be the disjoint union of the Gx-, x E H . 

• For every edge (x, y) G H, add 2k edges between W^ and Wy and similarly 2k edges 
between W~ and W~ . This can be done deterministically in such a way that the 
resulting graph, which we denote H^ , is d-regular. 
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We write a spin configuration on H'^ or H'^ as a = {g^x)x<^H where o^ is the restriction of a to 
Gx- We write Y^ = Y{a^) for the phase of each a^, and 3^(a) = iy{o^))xeH ^ {0, 1}^. Write 
ZjjG (y) for the partition function for the two-spin model on H'-' restricted to configurations 
of phase 3^, and define hkewise Zfjciy)- 
Recalhng (4.1), let 

r = z{h^ ®^ h+)z{h' ®^ h'), e = z{h^ ®^ h-y, 

and note that for anti- ferromagnetic two-spin models in non- uniqueness regimes, B > F. 

Lemma 4.3. For G satisfying properties (l)-(lll), 

\(1 - e)/2]"' < ^HG/ZfjG < fl . A-rn 

LI J- Cj/^J — Y2k\E(H)\(QlYYkMA^-CVT{H) — ^ ^ ^> ■ 



Proof. By (ll), 

for all y G {0, 1}-^. By (ih), the ratio 



1 _ g)"' < 2"^ """^^^ < (1 + e)"^ (4.5) 

Zqg 



^^ ^e-f^2w^ {i,j)&E(HG)\E{HG) 

is within a (1 ± e)™ factor of 

Y. n ^^^ (^H/+) n ^^^^' ^^■)' 

which by direct calculation equals 

p2fc|_B(//)| /Q /■pN2fccut(3;) 

where cut(3^) = \{{x,y) G E{H) : Y^ ^ Yy}\, the number of edges crossing the cut of H 
induced by y. Combining with (4.5) gives 

^HO = E ^^.Z^aiy) < (1 + ef^T'''\Em(^Q/YfkMA.-ov.iH)z^^ 
and similarly 

Rearranging gives the stated result. D 

Using this lemma we now complete the reduction to approximate max-cut: 

Proof of Thms. 1 and 2. Let if be a 3-regular graph on m vertices, and note that the max- 
imum cut of H is at least 3m/4, the expected value of a random cut. Construct if*^, if^ 
as above. Since H^ is a disjoint collection of constant-size graphs, its partition function can 
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be computed in polynomial time. Suppose Z^a could be approximated within a factor of 
gc|// I jj-^ polynomial time for any c > 0: rearranging the result of Lem. 4.3 gives 

/ Zhg/Zqg \ , / Zhg/Zqg 

log T.-yi.\K(mw. , .N^ log 



^ ^ ' ^ < MAX-CUT(if) < ^^^ , , ^r.,T^^ ^ (4-6) 



2fciog(e/r) - ' '- 2A;iog(e/r) 

so within polynomial time one obtains upper and lower bounds for MAX-CUT(if ) which differ 
by ©[(cIGI + l)m/k]. Taking k large and c small then allows to compute MAX-cUT(if) up 
to an arbitrarily small multiplicative error: that is, we have completed the reduction to a 
PRAS for MAX-CUT on 3-regular graphs, in contradiction of the result of [1]. D 
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